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' Abstract. Previous work on the angular momentum balance and meridional 

\ circulation of the solar convective zone (SCZ) generally consists either of semi- 

, analytic approaches in which a simple turbulence model is adopted, or full direct 

[t I ' numerical simulation (DNS) of the hydrodynamics or magnetohydrodynamics 

(MHD). In both instances the inclusion of magnetic fields has been troublesome. 
Also, both approaches have had difficulty reproducing the known angular velocity 
profile of the SCZ; this is the Taylor Number Puzzle. 

I discuss preliminary work in which I incorporate magnetic fields into a 
^ ■ viscoelastic turbulence model for the SCZ and apply this to the problem of 

' meridional circulation and angular momentum balance. I suggest that such an 

\ approach may help solve the Taylor Number Puzzle of the SCZ and bring theo- 

i retical predictions for the large-scale motion of the SCZ in line with observations. 

(N 
O . 

\^ , The Taylor number puzzle is this: in mean- field simulations of the combined 

■ system of meridional circulation and rotation of the solar convective zone (SCZ), 

r~| . one may come reasonably close to reproducing the observed rotation profile by 

Qh' making the Taylor number artificially small. If on the other hand one uses 

Q . realistic values for the Taylor number, the meridional circulation is much larger 

' in magnitude than otherwise. This results in a theoretical angular momentum 

c/3 . profile that differs markedly from the observed profile. (Brandenburg et al. 199^) 

^ ' To be fair, the puzzle may have already been solved, at least to a degree. 

That is because the situation changes if one allows for various mechanisms that 
create a significant baroclinic vector, as baroclinicity appears as a source term 
in the meridional circulation. This baroclinicity could be generated by latitu- 
\ dinal variations in the thermal transport, anisotropy of the thermal diffusivity 

tensor, subadiabaticity of the solar t ac hocline, or some combinati on of these 
effects (jKitchatinov fc Riidigedll995al lbl: lR^mDe]ll2nn,ll : fMeisch et ahl EnOfi). An- 



other possibility is that additional Reynolds stresses due to the anisotropic ki- 
netic alpha (AKA) effect may b e important to getting theory to match observa- 
tion ( Rekowski k, Rudigeilll99^ ). Still, I take an alternative point of view that 



the solution lies in part neither in baroclinicity nor in Reynolds stresses, but in 
the physics of Maxwell stresses. Here I focus particularly on turbulent Maxwell 
stresses rather than integral-scale fields. 

The provenance of the Taylor number is laboratory fluid dynamics. The 
definition varies, but the most common expression given is that 

Ta=^^. (1) 



Ta is a dimensionless quantity used to characterize the relative strength of cen- 
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trifugal driving and viscous dissipation of meridional circulation, particularly in 
cylindrical and spherical Couette flow. This latter problem, the flow between two 
rotating spheres, bears some obvious similarities to the SCZ. In the laboratory, a 
meridional circulation is driven when the spheres are differentially rotating. The 
magnitude of this secondary flow relative to the primary azimuthal flow increases 
as Ta (or equivalently Re) increases. 

To understand what drives the circulation in both cases, let us look at the 
vorticity equation for a compressible fluid 



Actually, that is a bit of a lie; the viscous term for a compressible fluid is more 
troublesome than indicated here, which is why I feel compelled to use a ~ in- 
stead of = above. Let us not worry about these details for the moment. The 
viscous term damps meridional circulation^. In the case of the laboratory flow 
one could argue that viscosity also drives circulation through viscous coupling 
to the spheres, creating an adverse angular momentum profile; that is not our 
concern here. In the SCZ there are also contributions to the Reynolds stresses 
that are often referred to as non-diffusive viscous terms (this terminology reflects 
a broader sense of the word "viscosity" than I am adopting here). These stress 
components are very important to the inainte nance of differential rotation in the 
SCZ; this is the lambda-effect (iRiidigeJIlOSgT ). The primary stresses due to the 
lambda-effect are r(j) and 9(j) stresses, which do not appear in the (p component of 
the vorticity equation. Rather, the influence of the lambda-effect on meridional 
circulation is indirect. The lambda-effect bears on the meridional transport of 
angular momentum, which contributes to the nonuniform angular velocity pro- 
file, which in turn generates a meridional circulation through the V x (n x a;) 
term. This is clearly important for the SCZ, but again, it is not our concern 
here. 

The physical nature of the advection/stretching term, the first term on the 
right-hand side (RHS) of eq. ©, is best seen in cylindrical coordinates. Then, 
in an inertial frame of reference, when the meridional circulation is negligible. 



and this in turn is simply the curl of the specific centrifugal force in the locally 
corotating reference frame. This tells us that there is a source of meridional 
circulation if the rotation is not constant on cylinders. 

For the measured rotation of the SCZ, dz^{R,z) <^ djift{R, z); the RHS of 
eq. ((31) is dynamically significant. Under such circumstances, a barotropic fluid 
cannot undergo pure rotation and be in simultaneous gravitational, centrifugal, 
and pressure equilibrium. There must be a meridional circulation. The sense of 
action of the centrifugal driving can be read directly off of graph s of isostrophic 
contou rs. By inspection of the results of helioseismology, e.g. ISchou Sz et al.l 
(^9^, the predominant sign of 9^(r2^) in the northern hemisphere is negative. 



dtu V X {u X uj) + — V/3 X VP + uV'^uj 





(3) 



^Strictly speaking, this depends upon how broadly one defines viscosity; here I am using the 
term in a narrow sense. 
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The centrifugal force should drive a counterclockwise (CCW) circulation in the 
usual notation (that is, towards the poles at the surface and towards the equator 
at the base of the SCZ). This is similar to the driving of a secondary flow in 
spherical Couette flow when the inner sphere is rotating faster than the outer 
sphere; at moderate Re this flow is dominated by a single large axisymmetric cell 
in each hemisphere, CCW in the northern hemisphere. 

Applied to the SCZ using the turbu l ent vi scosity for i^, the numerical value 
of Ta is very high; iRekowski fc Eiidiged (|l998l 'l suggest Ta ~ 10^^^. Turbulent 
viscosity is simply not effective at dampening meridional circulation. Absent 
additional physics in the SCZ, the angular velocity profile as measured should 
strongly drive a meridional circulation. This circulation in turn should redis- 
tribute angular momentum, again through the V x (li x w) term, which in the 
SCZ should drive the system to a Taylor-Proudman state. This is not seen. 

We come now to the second term on the RHS in eq. ©. This term is of 
course irrelevant for the laboratory case of spherical Couette flow, but probably 
important for the SCZ. In principle, the baroclinic vector could balance out the 
centrifugal term in the meridional vorticity equation, creating a meridional flow 
much slower in magnitude than one would otherwise expect. Of course, one is not 
free to adopt whatever density and pressure profiles one wishes. In addition to the 
sought mechanical equilibrium, the system must also satisfy thermal equilibrium. 
Sadly, the direction and magnitude of energy flux is more difficult to ascertain 
in the convective zone than it is in the radiative zone. Simply note here that 
for the baroclinic vector to have the right sign to balance the centrifugal forces, 
following an isopycnic surface, pressure and hence entropy must increase as one 
moves from the equator to the poles. The temperature gradient must be less 
inclined to the axis of rotation than either the pressure or density gradients. 
One might expect this to cause an observable latitudinal temperature difference. 
Various means by which Nature may create such baroclinicity and thus solve the 
puzzle were briefly referenced at the start of this paper. 

As a hint for what is to come, note that the physics of turbulence is not 
included in eq. or at least not included very well. Granted, it is common 
to think of turbulence as creating an effective viscosity, a notion that I have 
already appealed to in writing the Taylor number for the SCZ in the first place. 
It is better, though, to write the effect of turbulence in eq. ^ applied to the 
SCZ as t he divergence of the tur bulent stress tensor, which includes Reynolds 
stresses. iKitchatinov et ah I (|1994') suggest an effect in convective turbulence in 
rotating, stratified fluids; this AKA-effect was mentioned previously. It creates 
rO and 96 Reynolds stresse s, and so is a source term for meridional circulation. 
iRekowski fc Riidiged (119981 ) show that this effect can help solve the Taylor num- 
ber puzzle as well. 

That is where things stand in the literature, in an extremely small nutshell. 
What physics might be missing, if any? Certainly one piece of missing physics is 
that the sun can and does support magnetic fields. Magnetic fields in turn create 
a stress, the Maxwell stress, the divergence of which is a force that acts upon 
the fluid. None of this is taken into account in eq. 0. Naturally, it has long 
been recognized that magnetic fields might be dynamically significant to the solar 
interior. The problem is that magnetic fields are difficult to treat. Here I am 
going to discuss some recent work on a way of treating one particular component 
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of the magnetic field — the anisotropic, smah-scale, turbulent magnetic field — 
and I am going to explain what I think is a very interesting reason to suspect 
that this might help get the angular velocity profile of the SCZ right. 

Let me take what appears to be a detour by returning to spherical Couette 
flow. Even an incompressible fluid in spherical Couette flow, for the case we have 
considered where the inner sphere rotates faster than the outer, can produce a 
reversed, clockwise (CW) flow rather than a CCW flow, without the help of any 
baroclinicity. All that is required is a viscoelastic flu id rather than a viscous one 
( Yamaguchi et al1ll997l : lYamaguchi &; MatsuilllQQ?! ^. In such a fluid, there is an 
additional source of meridional vorticity that is not captured in eq. HI 

A viscoelastic fluid can be thought of as a viscous fluid where the stress has 
a decidedly finite (non-zero) relaxation time s. On short timescales, the fluid 
behaves elastically; on long timescales the fluid behaves viscously. Such behavior 
is typically a result of the stress being due to long polymers rather than colliding 
point-like particles. For example, bread dough made from wheat is viscoelastic 
because of the gluten polymer; likewise, the ovalbumin protein (a polymer) makes 
egg whites viscoelastic. In the case of steady shear, a viscoelastic fluid produces 
a streamwise elastic stress, creating a positive (first) normal stress difference. ^ 

The primary (azimuthal) flow in spherical Couette flow is a type of of circu- 
lar viscometric ^ flow. Almost universally, viscoelastic fluids show this remark- 
able property in rotating viscometric flows: when viscoelastic effects dominate, 
there is a general tendency for the secondary (meridional) flows in such fluids 
to be opposite to the secondary flows that would have been obtained for the 
centrifugally-driven circulation of a purely viscous fluid. This is known as sec- 
ondary flow reversal. It should be clear that a similar effect in the SCZ would 
slow the inertially-driven meridional circulation, which is what we need. 

For this picture of secondary flow reversal to be applicable to the SCZ, 
we require the dynamics of turbulence in the SCZ to be similar to the dynam- 
ics of a viscoelastic fluid. But that is no problem. Both kinematically and 
dynamically, a turbulent, tangled magnetic field bears more semblance to a net- 
work of polyn iers than it does to co lliding p oint-l ike particles flOeilvidl l200ll : 
IWilliamsl l2nnil : iLongcoue etini2nn.'j : lOgilvid 1200^ : IWilhamsl 120041 '! . Let over- 
bars denote averages and primes denote fluctuations. Neglecting the technical 
issues of Reynolds versus Favre averaging, the turbulent stress Wij can be split 
into a Reynolds stress Rij and a turbulent Maxwell stress Mij, this latter being 
the sum A4ij = Mij — \Mkk^ij where AirMij = B'-B'j. Let us suppose that the 
Reynolds stress will consist of a turbulent viscosity and lambda-effect terms, as 
in previous mean-field work, and suppose as well that the turbulent magnetic 
field will contribute viscoelastic behavior, like a network of polymers. In princi- 
ple, pressure terms from both types of stress may contribute as well, but I am 
ignoring isotropic turbulent stresses here. 



^The first normal stress difference A*'i is the difference between tfie streamwise normal stress and 
tlie cross-stream normal stress and is zero for a Newtonian fluid. 

■^A viscometric flow is a simple shear flow so that the relative-right Cauchy- Green tensor can 
be written as the sum of the unit tensor and the first and second Rivlin-Eriksen tensors. An 
axisymmetric fiow with negligible meridional circulation is one example. 
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It is the normal stress difference A'^i (due to tlie elastic streamwise tension 
stress) of a viscoelastic fluid under steady shear that is largely responsible for 
the secondar y flow reversal in viscometric flows such as our laboratory example. 
According to IWilliami ()2004l ^ , such a streamwise tension stress should also exist 
in the SCZ, corresponding to the statistical alignment of the tangled field with 
the (p direction. Note that for the purposes of determining the stress, there is 
no distinction between the vectors B and —B. So, it is quite possible to have a 
mean stress, even an anisotropic one such as I am considering here, without a 
mean field. 

Formally, the streamwise stress in viscoelastic models for the turbulent mag- 
netic stress comes through a term proportional to the second Rivlin-Eriksen 
tensor. This tensor is quite simple, all components being zero except for the 
azimuthal normal component, i.e. the component. In a laboratory fluid, this 
stress, and so the normal stress difference A'^i as well, is equal to the product 
of viscous stress i^y and the Weissenberg number 57, although strictly this de- 
pends somewhat on how one chooses to define the relaxation time s. Here the 
situation is slightly different. Physically it makes sense to treat the turbulent 
magnetic field and the turbulent velocity fluctuations on different footings, and 
this complicates matters. There is now a total viscous stress which is the sum 
of the Reynolds viscous stress and the magnetic viscous stress, and there is an 
overall effective relaxation rate which combines the Reynolds stress relaxation 
rate and the turbulent Maxwell stress relaxation rate. In order to keep things 
simple here, let us just say that if we are not too far away from equipartition, 
the effect of viscoelasticity is to introduce a stress that is equal to usj'^ times 
some constant of order unity. 

The curl of the divergence for such a stress is potentially nonzero and is a 
source for the meridional circulation. For a stress of the form of an axisymmetric 
dyad acpcp such as we have here, the curl of the divergence is most simply expressed 
in cylindrical coordinates: 

V X (V • {a4>4>)) = -d, (^^a) (4) 

In the northern hemisphere in the mid-latitudes, the tendency is for the overall 
magnitude of the shear to decrease with z. Since the elastic stress is a positive, 
tension stress, the result is that the turbulent elasticity is a source for positive 
vorticity, which then drives a clockwise circulation. This is the correct direction 
to balance the inertially-driven CCW circulation, in principle. 

Let us estimate the stress needed to perform this balancing act and see if it 
is reasonable. One way of describing the magnitude of the required stress is to 
compare it to the viscous stress. Using a reference value of = 10^^'^ cm^/s and 
7 ~ Rn' ~ R{An)/{AR) ~ IQ-'^-^/s, the viscous stress is -yu ~ 10^ cm^/s^. The 
ratio of the required turbulent streamwise stress to the turbulent viscous stress 
is then 

^ = .7-10^. (5) 
z^7 

A current estimate for the ratio of streamwise stress to viscous stress for turbu- 
lence driven by the magnetorotational instability (MRI) in accretion di sks is that 
it is of the order of lO'''^ to 10^, based on shearing-sheet simulations ( Williani^ 
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The interior of the SCZ is another matter, but I suggest that it is encour- 
aging that the estimate in eq. © differs by no more than one and a half orders 
of magnitude from the MRI results. 

Alternatively, the strength of the turbulent, toroidal magnetic field can be 
expressed in terms of the Alfven velocity, v\ = B'^ /{Airp). We want the centrifu- 
gal source of negative (CCW) vorticity, Rdz{^^), to be balanced by the elastic 
Faraday tension source of positive (CW) vorticity, —R^^dz{v\). Consider a point 
of mid-depth in the convection zone, somewhere in the mid-latitudes. Here 

Rdzin^) Rn^ (6) 

and 

Using the fiducial values of = 27r x 400 nHz, AO = 27r x 20nHz, R = Rq, the 
corresponding azimuthal Alfven velocity is va — few x lO^cm/s, or of the order 
of f rot /few. For a characteristic density in the middle of the convection zone 
oi p = 0.05 g/cm^ this results in fields of roughly 10^'^ kG. This does not seem 
unreasonable. 

I conclude, then, that simulations of the SCZ that include the dynamic 
effect of the Faraday tension of a toroidal field — such as through a viscoelastic 
turbulence model, direct numerical simulation of the field, or otherwise — may 
find that this makes a significant difference in the meridional circulation and 
angular momentum balance; in particular, I hypothesize that the presence of 
a magnetic field may make the calculated angular velocity profile more nearly 
coincide with the observed profile, and be an important ingredient in the solution 
of the Taylor number puzzle. 
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